%[y, Fs] = audioread("1皮鞋走路.mp3")
% y 音频数据矩阵（列数为声道数，如单声道：[N, 1]，立体声：[N, 2]）
% Fs 采样率（Hz，如 44100）

clc
clear -regexp ^(?!Num$).%去除除了Num之外的数据，Num用于滤波
Fs=44100;
[a1,b1]=text_1(1);%a1是原音频，b1是噪音音频,c1是去噪音频
c1=filter(Num,1,b1);%获得去噪音频数据矩阵
audiowrite('c1.mp3', c1, Fs);

[a2,b2]=text_1(2);%a1是原音频，b1是噪音音频,c1是去噪音频
c2=filter(Num,1,b2);%获得去噪音频数据矩阵
audiowrite('c1.mp3', c2, Fs);

[a3,b3]=text_1(3);%a1是原音频，b1是噪音音频,c1是去噪音频
c3=filter(Num,1,b3);%获得去噪音频数据矩阵
audiowrite('c1.mp3', c3, Fs);



figure;
%****************************
% 音频1的时域图
subplot(3, 2, 1);
t1 = (0:length(b1)-1) / Fs;
plot(t1, b1);
xlabel('时间');
ylabel('振幅');
title("音频1");
grid on;
%音频1去噪的时域图
subplot(3, 2, 2);
t1_1= (0:length(c1)-1) / Fs;
plot(t1_1, c1);
xlabel('时间');
ylabel('振幅');
title("音频1去噪");
grid on;

%****************************
% 音频2的时域图
subplot(3, 2, 3);
t2 = (0:length(b2)-1) / Fs;
plot(t2, b2);
xlabel('时间');
ylabel('振幅');
title("音频2");
grid on;
%音频2去噪的时域图
subplot(3, 2, 4);
t2_1 = (0:length(c2)-1) / Fs;
plot(t2_1, c2);
xlabel('时间');
ylabel('振幅');
title("音频2去噪");
grid on;

%****************************
% 音频3的时域图
subplot(3, 2, 5);
t3 = (0:length(b3)-1) / Fs;
plot(t3, b3);
xlabel('时间');
ylabel('振幅');
title("音频2");
grid on;
%音频3去噪的时域图
subplot(3, 2, 6);
t3_1 = (0:length(c3)-1) / Fs;
plot(t3_1, c3);
xlabel('时间');
ylabel('振幅');
title("音频2去噪");
grid on;
%时域图对比结束********************************


%***********************************************


%幅频图***************************************
%音频1的幅频图*********************************
figure
subplot(3, 2, 1);
%对信号序列进行dft变换并取绝对值
N = length(b1);
Y = fft(b1, N);
Y_abs = abs(Y); 
%将横坐标转换为频率（hz）
f = (0:N - 1) * Fs / N;
%绘制幅频响应图
plot(f, Y_abs);
xlabel('频率');
ylabel('幅值');
title('音频1幅频图');
%音频1去噪******
subplot(3, 2, 2);
%对信号序列进行dft变换并取绝对值
N = length(c1);
Y = fft(c1, N);
Y_abs = abs(Y); 
%将横坐标转换为频率（hz）
f = (0:N - 1) * Fs / N;
%绘制幅频响应图
plot(f, Y_abs);
xlabel('频率');
ylabel('幅值');
title('音频1去噪');
%****************************************

%音频2的幅频图****************************************
subplot(3, 2, 3);
%对信号序列进行dft变换并取绝对值
N = length(b2);
Y = fft(b2, N);
Y_abs = abs(Y); 
%将横坐标转换为频率（hz）
f = (0:N - 1) * Fs / N;
%绘制幅频响应图
plot(f, Y_abs);
xlabel('频率');
ylabel('幅值');
title('音频2幅频图');
%音频2去噪******
subplot(3, 2, 4);
%对信号序列进行dft变换并取绝对值
N = length(c2);
Y = fft(c2, N);
Y_abs = abs(Y); 
%将横坐标转换为频率（hz）
f = (0:N - 1) * Fs / N;
%绘制幅频响应图
plot(f, Y_abs);
xlabel('频率');
ylabel('幅值');
title('音频2去噪');
%****************************************

%音频3的幅频图****************************************
subplot(3, 2, 5);
%对信号序列进行dft变换并取绝对值
N = length(b3);
Y = fft(b3, N);
Y_abs = abs(Y); 
%将横坐标转换为频率（hz）
f = (0:N - 1) * Fs / N;
%绘制幅频响应图
plot(f, Y_abs);
xlabel('频率');
ylabel('幅值');
title('音频3幅频图');
%音频2去噪******
subplot(3, 2, 6);
%对信号序列进行dft变换并取绝对值
N = length(c3);
Y = fft(c3, N);
Y_abs = abs(Y); 
%将横坐标转换为频率（hz）
f = (0:N - 1) * Fs / N;
%绘制幅频响应图
plot(f, Y_abs);
xlabel('频率');
ylabel('幅值');
title('音频3去噪');
%****************************************